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Abstract 

We study the localization of a random heteropolymer onto an homogeneous 
surface, the problem which is equivalent to the wetting of an interface at disor- 
dered substrate in two dimensions, via replica trick by using the Green's func- 
tion technique. The exact treatment of one- and two-replica binding states is 
used to compute the free energy of the random heteropolymer. We present 
analytical results for two particular cases: (i) nearly statistically symmetric 
copolymer in the vicinity of the threshold of the annealed problem, and (it) the 
asymmetric polymer with the interaction part of the annealed Hamiltonian 
being nearly zero. In both cases the localization is due to two-replica binding 
states. In the case (i) the two-replica binding state exists both above and 
below of the one-replica binding state. In the case (ii) the energy of the two- 
replica binding state at the transition is finite. A schematic phase diagram 
of the localization-delocalization transition of the random heteropolymer is 
suggested. 

PACS numbers: 05.20.-y, 61.41. +e, 64.70.-p, 72.15.Rn 
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Adsorption of a polymer chain onto a surface is of large practical interest ranging from 
biological physics to technological applications JTJ- (and citations therein). The adsorp- 
tion of a polymer with heterogeneous sequence structure is relevant in connection with the 
study of the behavior of proteins near surfaces. The problem of wetting in two dimensions 
fT0|j - |fn}| is closely related to the adsorption of a symmetric heterogeneous polymer onto 



a surface. The polymer adsorption has been studied in connection with denaturation of 
double-stranded DNA in solution [TjJ- [[HJ and growth problems [ |17| . 

The role of disorder on the absorption of the random heteropolymer remains a subject 
of controversy. This question has been addressed recently by several groups |TT[-^]. The 
quenched part of disorder is concluded to be irrelevant in |TT| , |T2"|j , while on the contrary, 
the work fl3| shows that the quenched part of disorder shifts the transition temperature 



of the localization. Very recently the similar problem was considered in WM by using a 



functional renormalization group. These authors predict that the localization transition for 
the symmetric heteropolymer is of Kosterlitz-Thouless type. 

In this Paper we study the adsorption of a random asymmetric heteropolymer via replica 
trick by using the Green's function technique. In contrast to the previous studies the Green's 
function method allows an exact consideration of one-replica (lp) and two-replica (2p) bind- 
ing states. Our analysis shows that for statistically symmetric heteropolymer the 2p binding 
state exists at the localization transition of the lp binding state. The exact solution of the 
two-replica problem at the point where the interaction part of the annealed Hamiltonian 
is zero, indicates that the localization transition is first-order. We incorporate the one- 
replica and two-replica localized states into a novel and heuristic procedure to compute the 
quenched free energy. 

The partition function of a polymer containing N segments interacting with the surface 

is 



Dz(s) exp 



- ./„ V ds 



1 ( N ds (^i \ _ / ,/.sv;,u(.s))- / ( i^(,)v llll ( z{ ,)) 



where Vq{z) = (oo, z ^ 0; 0, z > 0) is the repulsive interaction potential with the wall, I 
is the statistical segment length. The heterogeneity of the polymer is described by random 
Gaussian variables C( s ); which are characterized by the moments ((s) = (, and ((s)((s r ) = 
Q 2 + A5(s — s'). If C — 0, then the number of the monomers which are attracted or repelled 
from the well is on average the same, so that the heteropolymer is statistically symmetric. On 
the contrary, if ( ^ 0, there is an excess of the monomers, which are repelled from [Q > 0) 
or attracted to (£ < 0) the well. In this case the heteropolymer is asymmetric The 
attractive interaction with the surface will be modelled by the potential Vi nt (z) = u5(z — z°), 
where z° is small but nonzero (see below). This choice of the potential allows the exact 
treatment of the problem with two replicas, as will be shown below. 

Imaging s to be an axis perpendicular to z we interpret Eq.(|l]) as a partition function 
of a directed line interacting with a heterogeneous substrate at z — 0, which is the wetting 
problem. The random variables ((s) are now attributed to the substrate. The wetting 
interpretation of the polymer problem is a particular case of the relation between the polymer 
in d dimensions and directed polymer in d+1 dimensions. In the case of wetting the condition 
£ 7^ means that on average the interface interacts with the substrate. If £ > 0, then the 
interface is repelled from the substrate. 
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Performing the average over ((s) by using the replica trick we obtain the replica partition 
function as 



Z n = Z» = J Dz a (s) ex P (- J ds ^^T^ 2 + v oM*)) + <y**M*))) + 

a - r N 

- V] / dsV int (z a (s))V int (z b (s))). (2) 

2 a,b=l J ° 

The free energy has to be computed as follows 

I— (3) 

Due to the fact that @ contains only the one-fold integral over s, the partition function 
Z n can be interpreted as the probability amplitude of a quantum mechanical system of n 
particles associated with the Hamiltonian 

n n 

H n = J2(~ D ^l + V { z a ) - (35{z a -z ))- Au 2 ]T 5(z a - z°)5(z b - z°), (4) 

a=l a<b 

where D = l 2 /2, (3 = Am 2 /2<5 — C) with So = 1/5(0) being of the order of magnitude 
equal to the width of the potential well, and ( = (u. In the following we will speak of 
particles instead of replicas. The first term in @ is associated with the annealed average 
of the free energy. The annealed part of H n is equivalent to the localization problem of a 
quantum mechanical particle in the vicinity of the wall, and can be solved exactly. The 
binding problem for (|]) with n = 2 can also be solved exactly. The Hamiltonian (|4]) at 
given n can be interpreted as a Hamiltonian of n polymers interacting with the surface. The 
2nd term in the first sum describes the attraction of the monomers to the surface, which is 
independent of each other. The second sum in (|j) gives an additional attractive interaction, 
if the monomers belonging to different polymers contact the surface simultaneously. We are 
not aware, if such an interaction can be realized in reality. Eq.(^j) at n = 2 and for (3 = and 
^o( z a) — is exactly the Hamiltonian of a quantum mechanical particle in two dimensional 
delta potential. As it is well-known from text books the binding state in this case exists for 
infinitesimally weak potential. While YlZ=i Za ) ^ oes no ^ P ossess the radial symmetry, 
the Hamiltonian (f|) at n = 2 and Vq( z a ) ^ does not correspond to a quantum mechanical 
problem in a radial symmetric two dimensional potential well with an impenetrable core 
at the origin. The Hamiltonian H n at n = 2, Vq( z a ) = and j3 = is related to the 



Poland-Sheraga model |L5|. The essential difference to the Poland-Sheraga model consists 
in the fact that all contacts occur at z = z , i.e. the model (||) neglects the wiggling of the 
zipped polymer pair. It is evident that the last term in (^) favors the localization of replicas. 

To this end it is convenient to consider the one-replica Green's function, G(zi, N; z®) = 
(5(z(N) — Zi)8(z(0) — Zy)), associated with the annealed part of the Hamiltonian (£|). The 
Laplace transform with respect to N (the variable iV plays the role of the imaginary time for 
the quantum mechanical particle) of the perturbation expansion of G(zi, N; zf) in powers 
of the attraction strength (3 is a geometric series, which is summed as 
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rii o\ ry i o\ , a^oi Zi,p; z )G (z ,p; z^) , 

1 - (3G (z ,p; zq) 

where G (zi,p; z\) = (exp(— \z\ — z®\ ^/p/D) — exp(— \zi + z®\ ^/p/D)) /^ADp is the 
Laplace transform of the Green's function of the diffusion equation in the half space (z ^ 0) 
with the Dirichlet boundary condition at z — 0. The equation 

1 - PG (z ,p; z ) = l- (3{ADp)- l ' 2 {l - exp(-2z ^p~/D)) = 0, (6) 

which is the denominator of the 2nd term in the r.h.s. of Eq.(|5]), is the energy eigenvalue 
condition for lp (one-replica or one-particle) localized state. Identifying D as h 2 /2m, Eq.(||) 
coincides exactly with the eigenenergy condition for the localization of a quantum mechanical 
particle in an attractive Delta-potential placed at the distance z from the wall. The localized 
state corresponds to the solution of (|E|) p c > 0. The energy of the localized state is given 
by Ei t0 = —p c . It is easy to see from (|6]) that the localized state exists for (3 > [3 C = D/zq. 
The inverse Laplace transform of G(z ,p; z ) for weak binding (2z ^p/D ^ 1) and f3 > f3 c 
is obtained from (||) as 



G(z , N; z ) * ^(Vp c eMN Pc ) + ^= - exp(iV> c ) - ^ -J^z + 

^fp~ c exp(Np c ) eii(^/N~p c ) - exp(iVp c ) ert{y/Np c )), (7) 



D 

where p c = (/3//3 c — l) 2 D 3 / \f3 2 ] Zq) . In computing (^) we neglected the term proportional to 
S(N) i.e. we consider N in (|7]) to be positive. Taking into account the delta function in 
computing the two- replica binding state in Eqs. fl9[-|10|) leads to non significant changes. 

To study the effect of the non-diagonal part of H n in the case n = 2 we will consider 
the connected part of the two-replica Green's function G 2>c (zi, z 2 , N; z®, z 2 , 0) = (S(zi — 
z a (N))5(z 2 — Zb^N^S^z® — z a (0))5(z2— z 2 (0))) c , where a and b denote the replica indices (a ^ 
b). The perturbation expansion of the Green's function G^ c {z\, z 2 , N\ z®, z 2 , 0) in powers of 
the interaction (third and fourth terms in (P) is represented graphically in Fig.l. The dotted 
lines are associated with the two-replica interaction given in Eq.(^). The ends of the dotted 
lines are associated with z and the arc length Sj, which are ordered from the left to the right. 
An integration over Si has to be performed. Each part of the continuous line between two 
consecutive dotted lines is associated with the one-replica Green's function G(z°, s,; z°, Sj_i). 
The left (right) external lines are associated with G(z, N; z°, si) (G(z , Si; z°, 0)) (z be Zi, 
or Z2, while z° be z®, or z®). The graphical expansion in Fig.l, which visualizes the effect 
of two-replica interaction in H 2 in terms of space-time (N = —it) trajectories, shows that 
both trajectories contact the surface at the same time. Thus, the return probability to 
have two consecutive contacts is the square of that for one particle. This suggests that the 
localization of two particles interacting according to (|J) is closely related to the localization 
of one particle in two dimensions. The integral associated with a graph in Fig.l is a folding, 
so that the Laplace transform with respect to iV reduces the perturbation expansion in Fig.l 
to a geometrical series, which is summed as 

n i o on G 2 ( zi, z 2 ,p; z , z )G 2 (z , z ,p; z±, z 2 ) 

G 2 , c {zuz 2} p;z u l} z^) = a ~- , (8) 

1 - aG 2 {z ,z ,p;z ,z ) 
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where a = Au 2 , and G-iizxi z 1 ,P] z i, z 2 ) * s the Laplace transform of the product of two one- 
replica Green's functions G(zi, N; z1)G(z2, N; z®)- The denominator on the r.h.s. of Eq.(^) 
gives the eigenvalue condition for the two-replica (particle) bound state. We did not succeed 
to analyze the latter analytically, so that now we will consider the following particular 
cases: (i) approximately statistically symmetric polymer in the vicinity of the localization 
threshold of the annealed problem, and (ii) the asymmetric polymer under condition that 
the interaction part of the annealed Hamiltonian is approximately zero. The latter case can 
be realized by tuning the asymmetry parameter (. The two-replica partition function for 



directed polymers with random interactions was previously studied in [21]- | 22]j . 

To ensure the existence of Laplace transform of the one-replica Green's function squared, 
we introduce a short-time cutoff by replacing 1/N by 1/(N + a). The cutoff a along the 
polymer can be eliminated in favor of the transversal length ao via a = a^/4:Dh. The 
eigenvalue condition for the two-replica localized state, which consists in equality of the 
denominator of Eq.(§) to zero, is obtained for small p c , i.e. in the vicinity of the localization 
transition of the annealed problem, as 

= 1 - -^TM e M4p/iD)V(0,a 2 p/iD) + 

71 Zq p 

l>: ' l)<> //_2rr + ^exp(a^/4D)r(0,a^/4D) V (9) 



7TZ 2 (3 2 V. VD 

where T(0,x) is the incomplete gamma function. Eq.(|9|) yields for p in the vicinity of the 
one-replica binding transition, i.e. for small p c 

2 o2 



4£>exp(-7) tx z^ (3 



P2,c= T " exp(— -j^), (10) 



where 7 is the Euler number. The energy E 2)C of the two-replica bound state is —p2, c - For 
the symmetric case, a = 2/55 , P2,c decreases with decreasing So at fixed (3 ~ (3 C and z . This 
ensures the validity of the condition Iz^pj D ^ 1 (weak binding) that we used to derive 
Eq . (|I~0"D . A similar analysis slightly above the annealed threshold also results in Eq. fllO|) , so 
that the two-replica bound state exists both below and above the threshold of the one-replica 
bound state. The result (|10| ) shows that the two-replica bound state already exists at the 
one-replica localization transition. This is very reasonable and can be explained qualitatively 
as follows. The individual interactions with the surface contained in the annealed part of the 
Hamiltonian @ result in an increase of the probability to find the monomers of the polymer 
pair in the vicinity of the surface. This compensates the decrease of the probability, which 
is due to the wall potential J2a=i z a)i anc ^ thus shifts the threshold to the lower values. 

The asymmetric case (ii), (3 — > ; where no one- replica localized states exist can be 
realized by tuning the asymmetry parameter ( as it is seen from the definition of j3. For 
small j3 the lp Green's function may be approximated by its bare value, so that we obtain 



1 Notice that the case, when only the first term in ( JTT| ) is present, corresponds to localization of 
a QM particle in a shallow 2d potential well. The above procedure gives an exact solution of the 
problem, if one identifies the length oq with the width of the potential well. 
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(11) 



As above we replace 1/N in the first term on the right-hand side of flTT|) by 1/(N + a), where 
the cutoff a along the polymer can be eliminated in favor of the length a® via a = a^/AD. 
The zo-dependent terms in (|TT| ) are due to the boundary condition at the surface z = 0. In 
the case of adsorption onto an interface only the term l/(4irDN) will appear in Eq. (|TT|) , 
so that in this case the two-replica bound state will exactly coincide with that in a shallow 
potential well in two dimensions. Using (|TT|) we obtain from (|8|) the eigenvalue condition for 
the two-replica bound state as 

1 + 40^(4 VW^) - 2a x K Q {4 ^WR) ~ «i exp(p)r(0,p) = 0, (12) 

where K (x) is the modified Bessel function of the second kind, p = pa^/D, a± = a/inD, 
a = clo/zq. Identifying the cutoff a Q with z gives a = 1. It appears that the results are 
not sensitive to the choice of do- The numerical analysis of Eq.(p!2|) for different a (= 0.8, 1, 
1.2) yields the critical value a\ for the localization transition (= 0.5825, 0.7857, 1.085). It 
appears that at the transition pp 2c has the finite value (= 0.0001, 0.00195, 0.01963), i.e. the 
binding transition is a first-order transition. Above the transition, a\ > a\, there are two 
solutions for p 2c : f>2,c < P2 c an d Pi,c > P2 c - According to the ground state dominance argu- 
ment the larger value governs the behavior of the polymer for large N. The approximative 
consideration based on the Taylor expansion of the eigenvalue condition fll2"|) for small p 

l-ai(7 + ln2-21na) + a 1 p(-l+7 + 81n2/a 2 + lnp) + ... = (13) 

is in agreement with the results of numerical consideration of Eq. (|I2"D . The reason of the 
unusual first-order transition is due to the term phap in Eq.([H|). The latter is responsible 
for that the l.h.s. of fllTf ) has a minimum at finite p, which leads to the first-order transition. 
The physical reason of the first-order transition is due to the boundary condition at the wall 
z = 0, which results in a reduction of number of conformations, and drives the transition to 
be first-order. As we stressed above the problem associated with the Hamiltonian H n given 
by Eq.(|4]) for number of replicas n = 2 and j3 = is closely related to the Poland-Sheraga 
model [T5| , which was recently studied in JT8"| - PD| by taking into account the excluded- 



volume interaction between denaturated loops and the rest of the chain. The excluded 
volume interaction reduces the number of conformations and drives the transition to the 
first order. This is similar to the two-replica localized state, where the reduction of the 
number of conformations is due to the effect of the wall potential Vo(z). 

We now will consider the computation of the free energy by using the replica formula ([|) 
under taking into account the one-replica and two-replica binding states. In the case, if only 
lp (one-replica) binding states exist, Eq.(Q) gives straightforwardly — /3F = p c . However, 
the situation is nontrivial when a 2p (two-replica) binding state exists. Taking into account 
the two-replica states in the two-pair approximation we obtain the partition function Z n as 

Z n = Zl + n(n ~ 1] Z 2 , c Zr 2 + ... = exp(np c N) + e xp(p 2<c N + (n- 2)p c N) + 

(14) 
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where we have taken into account that Z\ and Zi, c behave for large N as exp(p c N) and 
exp(p2, c N), respectively, when both lp and 2p bound states exist. If no lp bound state 
exist, then Z\ = 1. The next terms in ( |T4| ) contain contribution of ternary and higher pairs. 
Using (||) and fllil) , the free energy is obtained to be proportional to exp((p 2 , c — 2p c )iV), 
hence the free energy is not extensive quantity for p2 )C — 2p c > . This shows that the 
two-pair approximation fllil) is insufficient for computing the free energy. The problem is 
due to the fact that the exponentially increasing term in (|TJ|) that originates from the two- 
particle bound state does not contain the factor n, as it is the case for the lp bound state. 
The factor n in the exponential of the latter ensures that it disappears in the limit n — > 0. 
To overcome the difficulty we suggest to take into account in the expansion of Z n the term 
containing the maximal number of unconnected pairs. Then, instead of ( |14|) we obtain 

Z 2n = Z l +2 - - Z 2c +..., (15) 

where we consider even n. The factor 2~ n T{2n + l)/r(n + 1) in Eq . fll5|) is the analytical 
continuation of (2n — 1)!! for arbitrary n. The free energy is now obtained from Eq.(^) as 

-J^-=Pc+^P2,c+- (16) 

The localization length £ ioc can be computed by using p 2;C as ^ oc ~ (p 2 , c /-D)~ 1//2 - The regime 
(ii) considered above is obtained from ([16]), if one puts p c = 0. Notice the plus sign in front 
of the 2p energy in ([16]). The computation of the 2p energy in the two-pair approximation 
(|TJ]) gives in the limit n — > instead of (|TB[) the minus sign in front of the second term. 
The inconsistency of this result is clearly demonstrated by the non extensiveness of the 
free energy computed by using Eq.(^). The procedure of taking into account the terms 
with maximal number of pairs can be justified by the following ground state dominance 
like argument: the term, exp(np2, c N), dominates over exp(p2,c-^) for large N and n ^ 1. 
The condition n ^ 1 is demanded in the procedure of introduction of the replica trick by 
considering the partition function Z n with n being positive integer. For the peculiarities of 
the limits iV — >• oo and n — > in the replica treatment of the directed polymer in disordered 
media see ||23|| . The extensiveness of the free energy is a posteriori justification of the above 
heuristic procedure. 

We now will use the results of the study of lp and 2p bound states of the replica Hamilto- 
nian (01) to construct the phase diagram of the localization of the heteropolymer in variables 
C and «i = Au 2 /47rD, which is shown schematically in Fig. 2. We know the behavior at a\ 
and af from the study of the cases (ii) and (i). The latter corresponds to the localization 
transition of the annealed Hamiltonian (lp bound state) in the symmetric case (( = 0). The 
dotted straight line is the localization line of the annealed Hamiltonian; £ c = —D/zq is the 
value of the asymmetry parameter in the limit a.\ — > 0. Since the random heteropolymer 
can arrange at the surface in such a way that pieces of the polymer which are attracted to 
the surface are in contact with the latter, while the pieces which are repelled are in the loops 
and tails (see 0] for a related discussion), the random quenched heteropolymer is expected 
to localize easier than the homopolymer or even easier than the annealed heteropolymer. 
Due to this the quenched part of the Hamiltonian (f|), which is responsible for 2p bound 
states, shifts the critical temperature to higher values, so that the localization line will cut 
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the horizontal axis on the left from the point af. The continuous line cannot end in the 
dotted line, while the analysis carried out in (i) in the vicinity of af applies also along the 
dotted line, so that the 2p localized state exists there. Thus, the continuous line will be 
on the left of the dotted line, and it will end in ( c . In order to make our prediction that 
the localization transition at a\ is first-order compatible with the prediction in (JTJJ that 
the transition at £ = is of Kosterlitz-Thouless type, one should assume that on the phase 
boundary in Fig. 2 there is a tricritical point separating the KT point at ( = with the 
first-order transition at a\. 

In summary, we have considered adsorption of a random heteropolymer onto a homoge- 
neous surface via replica trick by using the Green's function technique. We use the exact 
treatment of one- and two-replica binding states of the replica Hamiltonian to compute the 
free energy of the random heteropolymer. We have considered analytically two particular 
cases: (%) almost statistically symmetric polymer in the vicinity of the threshold of the an- 
nealed problem, and (ii) the asymmetric polymer where the interaction part of the annealed 
Hamiltonian is nearly zero. In the former case we have obtained that the localization of 
the polymer is due to the 2p binding state, which exists both above and below of the lp 
binding state, which corresponds to the localization transition of the annealed problem. We 
have obtained that in the case (ii) the energy of the 2p binding state at the transition is 
finite, i.e. the localization is a first-order transition. A schematic phase diagram of the 
localization-delocalization transition of the random heteropolymer is suggested. 
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Figure captions 

Fig.l The perturbation expansion of the connected part of the two-replica Green's func- 
tion. 

Fig. 2 The phase diagram of the localization-delocalization transition. 
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